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Abstract 

Using Korringa-Kohn-Rostoker coherent-potential approximation in the atomic- 
sphere approximation (KKR-ASA CPA) method for taking into account the effects 
of disorder, Gaspari-Gyorffy formahsm for calculating the electron-phonon coupling 
constant A, and Allen-Dynes equation for calculating Tc, we have studied the vari- 
ation of Tc in Mgi^xTaxB2 alloys as a function of Ta concentration. Our results 
show that the decreases with the addition of Ta for upto 40ai% and remains 
essentially zero from 60ot% to 80ot% of Ta. We also find TaB2 to be supercon- 
ducting, albeit at a lower temperature. Our analysis shows that the variation in T^ 
in Mgi^xTaxB2 is mostly dictated by the changes in the Bp density of states with 
the addition of Ta.. 



The experimental [1-4] and theoretical [5-12] efforts aimed at understanding 
the nature of superconductivity in MgB2 have made substantial progress since 
the discovery of superconductivity in it [1] . With an enhanced understanding 
of superconductivity in MgB2, attempts are now being made to understand 
the changes in the electronic structure and the superconducting properties 
of MgB2 upon alloying with various elements. Such efforts provide an op- 
portunity to explore the possibility of obtaining MgB2 alloys with improved 
superconducting properties. 



The changes in the superconducting properties of MgB2 due to substitution 
of various elements such as Be, Li, C, Al, Na, Zn, Zr and others have been 
studied experimentally [4,13-18]. Some theoretical work has also been reported 
within the rigid-band model. The main effects of alloying are seen to be (i) 
a decrease in transition temperature, Tc, with increasing concentration of the 
alloying elements except for Zn (and possibly Li) where Tc is seen to increase 
somewhat, (ii) changes in lattice parameters a and c, and (iii) a change in 
crystal structure. 
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With a view to understand the changes in the electronic structure and the su- 
perconducting properties as well as to clarify the reported superconductivity 
in TaB2 [19,20], we have carried out an ab initio study of Mgi_.j.TaxB2 alloys. 
We have used Korringa-Kohn-Rostoker coherent-potential approximation [21] 
within the atomic-sphere approximation (KKR-ASA CPA) method for taking 
into account the effects of disorder, Gaspari-Gyorffy formalism for calculating 
the electron-phonon coupling constant A, and Allen-Dynes equation for calcu- 
lating Tc in Mgi_xTaxB2 alloys as a function of Ta concentration. Such an 
attempt allows us to examine the possibility of superconductivity in TaB2, 
given the superconductivity in MgB2. We have analyzed our results in terms 
of the changes in the total density of states (DOS), in particular the changes 
in the B p contribution to the total DOS, as a function of Ta concentration. 

Based on our calculations, described below, we find that in Mgx-xTaxB2 alloys 
(i) the Tc decreases with the addition of Ta for upto 40 at%, remains essentially 
zero from 60at% to 80at%, and then rises to ~ 1.8 X for TaB2, (ii) the Tc 
for TaB2 is much lower than reported earlier [19], and (iii) the variation in 
is mostly dictated by the changes in the Bp densities of states as more and 
more Ta are added. Before we describe our results, we outline some of the 
computational details. 

The charge self-consistent electronic structure oi Mgi_xTaxB2 alloys as a func- 
tion of X has been calculated using the KKR-ASA CPA method [22] . We have 
used the CPA rather than a rigid-band model because CPA has been found to 
reliably describe the effects of disorder in metallic alloys. We parametrized the 
exchange-correlation potential as suggested by Perdew-Wang [23,24] within 
the generalized gradient approximation. The Brillouin zone (BZ) integration 
was carried out using 1215 k-points in the irreducible part of the BZ. For 
DOS calculations, we added a small imaginary component of 2 n)Ry to the 
energy and used 3887 k-points in the irreducible part of the BZ. The lattice 
constants for MgB2 and TaB2 were taken from experiments while for other 
compositions we used the Vegard's law. The Wigner- Seitz radii for Mg and 
Ta were shghtly larger than that of B. The sphere overlap, which is crucial in 
ASA, was less than 10% and the maximum I used was Imax — 3. 

The electron-phonon coupling constant A was calculated using Gaspari-Gyorffy 
[25] formalism with the charge self-consistent potentials of Mgi^xTaxB2 ob- 
tained with the KKR-ASA CPA method. Subsequently, the variation of T^. 
as a function of Ta concentration was calculated using Allen-Dynes equation 
[26]. The average values of phonon frequencies uin for TaB2 and MgB2 were 
taken from Refs. [27] and [28] respectively. For intermediate concentrations, 
we took uin to be the concentration- weighted average of MgB2 and TaB2. 
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Fig. 1. The calculated variation of superconducting transition temperature Tc with 
Ta concentration in Mgi^xTa^ B2, with n* = 0.09. 



The main results of our calculations are shown in Fig. 1, where we have plotted 
the variation in Tc of Mgi^,j.TaxB2 alloys as a function of Ta concentration. 
The calculations were carried out as described earlier with the same value of 
ji* — 0.09 for all the concentrations. The Tc for MgB2 is equal to ~ 30 K, which 
is consistent with the results of other works [28] with similar approximations. 
The corresponding A is equal to 0.73. As a function of Ta concentration we 
find that the decreases rapidly and goes to zero at around 40 af% of Ta as 
shown in Fig. 1. The Tc remains essentially zero until around 75 — 80at% of 
Ta, thereafter rising to 1.8 i^' for TaB2 with x — 1. The corresponding A for 
TaB2 is equal to 0.38.. 
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Fig. 2. The total density of states for Mgi^xTaxB2 as a function of Ta concentration 
X, calculated with the charge sclf-consistcnt KKR-ASA CPA method, as described 
in the text. Note that in (f), the vertical scale is different from the rest of the figures. 
The vertical line denotes the Fermi energy. 



In order to get further insight into the variation of Tc as a function of Ta 
concentration in Mg\-xTaxB2 alloys, we have analyzed the total DOS as well 
as the contribution of B p-electrons to the total DOS. In Figs. 2(a) — (/), we 
show the total DOS of A4gi_,j.TaxB2 alloys for Ta concentration ranging from 
a; = to a; = 1, calculated using the KKR-ASA CPA method as described 
earher. The DOS for MgB2 and TaB2, as shown in Figs. 2(a) and 2(/) are 
similar to previous calculations [28,29]. With small addition of Ta the Fermi 
energy, Ep, moves outward to accommodate the additional electrons resulting 
in a decrease in the total DOS at Ep for Mgi_xTaxB2 alloys, leading to 
a decrease in T^. Further addition of Ta allows Ep to be pinned in a valley 
around 50 — 60 at % of Ta, resulting in a decreased donor ability of the metaUic 
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a (a. u.) 


c/a 


N{Ef) 


B p 


MgB2 


5.834 


1.141 


3.331 


2.834 


Mgo.8Tao.2B2 


5.833 


1.124 


3.116 


1.785 


Mgo.&Tao.AB2 


5.831 


1.106 


0.466 


0.616 


Mgo.4Tao.QB2 


5.830 


1.088 


1.734 


0.663 


Mgo.2Tao.8B2 


5.829 


1.070 


2.791 


0.833 


TaB2 


5.828 


1.052 


4.980 


1.455 



Tabic 1 

The calculated lattice parameter a {a.u.), c/a, total DOS at the Fermi energy 
(states/Ry/atom) and B p partial DOS (states/Ry/spin). The lattice parameters 
for intermediate concentrations were determined using Vegard's law. 

plane {Mg — To) with the increase in the number of d electrons. Thus the 
overall contribution of B electronic states decreases as the concentration of 
Ta is increased up to around 80 atVo of Ta. Starting with around 80 at% of Ta 
the B contribution rises and so does the Tc as shown in Fig. 1. To substantiate 
these qualitative observations we have listed in Table I the total DOS and the 
Bp DOS at Ep. In addition, in Fig. 3, we show the Bp^ and pz partial DOS 
as a function of concentration. It is clear from Table I and Fig. 3 that the 
Bp DOS is responsible for the variation in Tc of Mgi_xTaxB2 alloys, and the 
possible loss of superconductivity in these alloys at around 60 at% of Ta can 
be attributed to a very small Bp^^y DOS at Ep. In Table I, we have also hsted 
the lattice constants used in our calculations. 



In conclusion, we have studied the variation of Tc in Mgi_xTaxB2 alloys as a 
function of Ta concentration. We have used the KKR-ASA CPA method for 
taking into account the effects of disorder and Gaspari-Georffy formalism for 
calculating the electron-phonon coupling constant A. The T^ is then calculated 
using the Allen-Dynes equation. Our results show that the T^. decreases with 
the addition of Ta for upto 40 at%, thereafter remains essentially zero up to 
80 at% of Ta. We find TaB2 to be superconducting, albeit at a low temper- 
ature. We have also shown that the variation in Tc is mostly dictated by the 
changes in the Bp density of states. 



We would like to thank Dr. I. A. Abrikosov for providing a copy of his locally 
self-consistent Green's function code. 
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Fig. 3. B p partial DOS showing (full line) and pz (dashed line) for Mgi^x Tax B2. 
The vertical line denotes the Fermi energy. Note that for x = 1 the vertical scale is 
different from the rest of the plots. 
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